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Abstract 

The physical properties of galactic winds are of paramount importance for our understanding of 
galaxy formation. Fortunately, they can be constrained using background quasars passing near 
star-forming galaxies (SFGs). From the 14 quasar-galaxy pairs in our VLT/SINFONI Mgll Pro¬ 
gram for Line Emitters (SIMPLE) sample, we reobserved the 10 brightest galaxies in Hct with the 
VLT/SINEONI with 0.7” seeing and the corresponding quasar with the VLT/UVES spectrograph. 
Applying geometrical arguments to these ten pairs, we find that four are likely probing galactic out¬ 
flows, three are likely probing extended gaseous disks, and the remaining three are not classifiable 
because they are viewed face-on. In this paper we present a detailed comparison between the line-of- 
sight kinematics and the host galaxy emission kinematics for the pairs suitable for wind studies. We 
find that the kinematic profile shapes (asymmetries) can be well reproduced by a purely geometrical 
wind model with a constant wind speed, except for one pair (towards J2357—2736) that has the small¬ 
est impact parameter 6 = 6 kpc and requires an accelerated wind flow. Globally, the outflow speeds 
are ^100 km s“^ and the mass ejection rates (or Mont) fhe gas traced by the low-ionization species 
are similar to the star formation rate (SER), meaning that the mass loading factor, 77 = Mout/SER, 
is ^1.0. The outflow speeds are also smaller than the local escape velocity, which implies that the 
outflows do not escape the galaxy halo and are likely to fall back into the interstellar medium. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: intergalactic medium — 
galaxies: kinematics and dynamics — quasars: absorption lines — quasars: indi¬ 
vidual (J0448+0950, J2357-2736, J0839+1112, J1441-h0443) 


1. INTRODUCTION 

Gurrently, under the Gold Dark Matter (GDM) sce¬ 
nario, galaxies form via the growth of initial den¬ 
sity fluctuations. This scenario is very successful be¬ 
cause the observed large scale structure is well matched 
by the cluster i ng of halos in N-body simulations (e.g. 
iSpringel et al.l I2QQ6F With the help of state-of- 
the art hydrodynani i cal simulations (iGenel et al.l |2Q14 
iMnratov et al.l 120151 : iSchave et al.l I2Q15F this scenario 
has recently successfully reproduced more complex ob¬ 
servables such as the galaxy m orphologies (|Genel et al.l 
|2Q14 IVogelsberger et alll2Ql4) . 

A major success of these recent hydro-simulations is a 
better understanding of disk formatio n at high-redshifts 
z > 1 . Indeed, these simulations (|Genel et al.l I 2 Q 12 L 
I2Q14D are in good agreement with the morphologies and 
kinematic observed for Lyman break galaxies, which ap¬ 
pear to be dominated by gas-rich rota ting disks, based 
on data from deep HST observations (|Elmegreen et al.l 
I2QQ7I: IWnvts et al.l 1201 ll) and large Integral Eield Units 
(lEU) surveys like SINS (lEorster Schreilper et al.|[2QQ6l: 
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iGenzel et al.l I 2 QQ 8 I: WE rster Schreiber et al.l l2QQ9bD and 


MASSIV (|Epinat et al.l 120121: iGontini et al.([ 2 Q 12 [) . 

One major problem remains, however, namely that 
the luminosity function for low mass galaxies {L < L^) 
is difficult t o reproduce . Eo r instance, the Illustris 
simulations (|Genel et al.l 120141) over-predicts the num¬ 
ber of 2 : = 0 low-mass galaxies d espite implementin g 
strong galactic outflows (but see ISchave et al.l I2Q15D . 
This tension can be recast in terms of the ‘galaxy for¬ 
mation efficiency’ (iMoster et al.ll2QlQl: iPapastergis et al.l 
l2Q12l:[Leauthaud et al.ll2QllLl2Q12l:lBehroozi et al.ll2QlQD . 

which is maximal for ^ galaxies and steeply de¬ 
creases in the low-mass regime {L < L^). In the 
low-mass regime, galactic winds, created by accumu¬ 
lated supernova explosions, are commonly invoked to 
expel baryons back into the inter-ga lactic medium 
(j White fc Erenkl[l991l: iDekel fc Silklll98^ since baryons 
i n these galaxy halos are expected to cool rapidly 
(IWhite fc ReesI 119781: IWhite fc Fren^ 119911: iKeres et all 
[20051: IDekel et al.ll2nn9D. 

Although galactic winds seem to occur in every 
star forming galaxy, their properties remain poorly 
constr ained despite many attempts at characterizing 
them ([Lehnert fc Heckmanl 119961: iHeckman et al.l I 200 QI: 


Martin I1998L 119991: iRupke et al.l 120051: iRnbin et alJ 

2 OIOI: Martin et al.l l2012f ). TTis lack of knowl¬ 


edge prevents us from correctly modeling galaxies 
in num erical simulations, which often require ad-hoc 
recip es (TOppenheimer fc Dayd120061: 1 Qppenheimer et al.l 
20101: iDubois fc Teyssierl 120081: iRoskar et al.l 2014 : 


Rosdahl et al.l 120131: iShen et al.ll2012L 120131: IDekel et aP 
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I2Q13D . In particular, the best estimates for the 
ejected mass rate (Mout) using standard galaxy ab- 
sorpti on lines are uncerta i n by o rders of magnitude 


(e.Qi. iHeckman et al. 

1199a 

2000 

; 

Pettini et al. 

2002; 

Martin et al.l 20021: 

Martini 

2005 

Martin et al. 

2012, 

20131). 


The main reason for the large uncertainties is that tra¬ 
ditional spectroscopy does not give information of the 
material physical location because the gas could be at a 
distance of 100 pc, 1 kpc, or 10 kpc from the galaxy. In¬ 
deed, the standard method usually uses the galaxy spec¬ 
trum and in some cases stacked galaxy spectra, to obtain 
the absorption lines corresponding to the outflowing ma¬ 
terials. However, background quasars have been recently 
used to constrain the prop erties of winds (|Bouche et al.l 
I 2 OI 2 I: iKacprzak et al.|[2014D using low-ionization absorp¬ 
tion lines, like Mgll (A2796, A2803). Indeed, when the 
quasar apparent location is close to the galaxy minor 
axis, the line of sight (LOS) is expected to intercept the 
wind. Thus background quasars give us the three main 
ingredients necessary for determining accurate ejected 
mass rates: the gas localization (impact parameter), the 
gas column density and the wind radial (de-projected) 
velocity, provided that the galaxy inclination is known. 

The background quasar technique also provides the 
ability to better constrain the ejected mass outflow rate 
and its relation to the galaxy star formation rate (SFR) 
via the so-called mass loading factor 77 = Mout/SFR, 
which is a critical ingredient for numerical simulations 
(iQppenheimer fc Davd I 2 QQ 6 I: iQppenheimer et all I 2 Q 1 QI: 
iDekel et ahl I2Q13D . As opposed to relying only on the 
galaxy spectra to study outflows, the background quasar 
method has several advantages: it gives us a more precise 
location of the absorbing gas relative to the galaxy and, 
because the quasar is seen as a point source, it also pro¬ 
vides us with a good characterization of the point spread 
function (PSF), an important ingredient for deriving the 
intrinsic galaxy properties from integral field unit (IFU) 
data. 

Recently, there has been progress in t his field with 
low-redshift z ^ 0.1 star-for ming galaxies (|Bouche et al.l 
I 2 QI 2 I: [Kacprzak et al.l 1201^ applying this technique. In 
this paper, we use the sample of 14 intermediate-redshift 
z ^ 1 galaxy-quasar pairs from the SINFONI Mgll Pro- 
gram for Line Emitters (SIMPLE, iBouche et al.l [200^ 
hereafter paper I) to constrain the outflow properties 
(e.g. mass ejection rate, outflow velocity) of star-forming 
galaxies when the quasar is suitably located relative to 
the foreground galaxy. 

The outline of this paper is as follows: § 2 describes 
the sample and the new VLT SINEONI/UVES data ac¬ 
quired. In § 3, we present the analysis of the SINEONI 
and UVES data together with the selection of pairs suit¬ 
able for wind studies (wind-pairs). In § 4, we describe 
our wind model and the derived outflow rates for the 
wind-pairs. We end with our conclusions and discussions 
in § 5. In this study we used the following cosmological 
parameters: Hq=70 km s“^, r2A=0.7 and f2M=0.3. 

2. THE SIMPLE SAMPLE 

Because the probability of finding galaxy-quasar pairs 
is very low, one must employ targeted strategies for gath¬ 
ering a suitable sample of galaxy-quasar pairs to study 


the properties of circumgalactic gas aro und galaxies 


which can lead to constraint s on outflows ( Bonche et al 


IKannrzak et a.lJ l2fnl) or inflows (|Bonche et al 
2013 : iPeronx et al.l 1201^ We thus designed the SIM¬ 


PLE survey to build a sample of intermediate redshift 

z 1 quasar-galaxy pai rs (paper I). _ 

The SIMPLE sample (|Bonche et al.ll2QQ7l) was selected 
with the following criteria: the rest-frame equivalent 
width of intervening Mg II absorptions detected 

in background quasar spectra had to be at least 2 A. 
This criterion ensures that the associated galaxies will 
be at small impact parameters (6 < 3” ), given the 
HA-impact paramet e r anti-correlation (e.g. iSteidell 19951: 
Bonch6 et al.l I2QQ6I: iMenard fc Chelonchel I2QQ9I: iChenI 


2Q12D . and thus that they will be located within the field- 


of-view of the lEU SINEONI (8” each side). Moreover, 
the absorber’s redshift must be 0.8 < 2 : < 1.0 so that the 
Ha emission line falls inside the SINEONI J band. These 
criteria led to the detection of 14 galaxies out of 21 (70% 
success rate) (|Bonche et al.l 1200"^ . 

The SINEONI data presented in iBonche et al.l (|2007[ ) 
were shallow with exposure times < 40 min and seeing 
conditions > 0.8”. Since we aim to precisely compare 
the host galaxy kinematics (derived from the Ha emis¬ 
sion line) with the kinematics of the absorbing material 
measured in the quasar line-of-sight, we acquired new 
VLT/SINEONI and VLT/UVES da^ . Erom the sample 
of 14 galaxies in IBouche et al.l (|2007[) . we re-observed a 
sub-sample of 10 galaxies, those with the highest initial 
Ha fluxes, with longer integration times (2-3 hr) and in 
better seeing conditions (< 0.8”). 

The SINEONI observations, done in service mode, 
were optimized by adopting a ‘on source’ dithering strat¬ 
egy designed to ensure a continuous integrat ion at the 
host location. The UVES (|Dekker et al.ll2Q0"Ql) data were 
taken in both visitor mode and service mode. 


2.1. SINFONI Data Reduction 

The d ata reduction was performed as in 
( 2007 b lEorster Schreiber et al.l (|2009a[) : 

( 2012), using the SINEONI pipeline (SPRED, 
iSchreiber et al.l l2004 I Abater et al.l I2006D comple¬ 
mented with cust om routines s uch as the OH sky line 
removal scheme of iDaviesI (|2007[ ) and the Laplacian e dge 
cosmic ray removal technique of Ivan DokkunJ (|2001f ). 

Regarding the wavelength calibration, we emphasize 
that we applied the heliocentric correction to the sky- 
subtracted frames, and each frame was associated with 
a single reference frame by cross-correlating each of the 
science frames spectrally against the reference frame (the 
first science exposure). Eor each observing block, we use 
the the quasar continuum to spatially register the vari¬ 
ous sets of observations. Einally, we created a co-added 
cube from all the individual sky-subtracted 600s expo¬ 
sures using a median clipping at 2.5cr. 

Elux calibration was performed on a night-by-night 
basis using the broadband magnitudes of the standards 
from 2MASS. The flux calibration is accurate to ^ 15%. 
Einally, the atmospheric transmission was calibrated out 
by dividing the science cubes by the integrated spectrum 
of the telluric standard. 

In Eigures [T] and [2] we present the flux, velocity and 
dispersion maps for each galaxy. 


Bonche et al. 

Bonche et al. 
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TABLE 1 

Summary of SINFONI 080.A-0364(B) observations. 


Field 

(1) 

-S^qso 

(2) 

(3) 

PSF(") 

(4) 

Texp (s) 

(5) 

Date 

(6) 

J0147+1258 

1.503 

4.025 

0.6 

9600 

2007-10-12 2008-01-03,04,09 

J0226-2857 

2.171 

4.515 

0.6 

9600 

2007-10-06 2008-01-03,05 

J0302-3216 

0.898 

2.27 

0.7 

7200 

2007-10-02 

J0448+0950 

2.115 

3.169 

0.8 

4800 

2007-12-04,16 

J0822+2243 

1.620 

2.749 

0.8 

4800 

2007-12-18 2008-01-03 

J0839+1112 

2.696 

2.316 

0.8 

4800 

2007-12-14,15,23 2008-01-01 

J0943+1034 

1.239 

3.525 

0.6 

7200 

2007-12-22 2008-01-06 

J1422-0001 

1.083 

3.185 

0.7 

9600 

2008-02-15,25 2008-03-14 

J1441+0443 

1.112 

2.223 

0.6 

12000 

2008-03-14,15,25 

J2357-2736 

1.732 

1.940 

0.6 

7200 

2007-10-02 


(1) Quasar name; (2) Quasar emission redshift; (3) Mgll rest-equivalent width; (4) FWHM of the seeing PSF; (5) Exposure time; (6) 

Dates of observations. 


2.2. UVES Observations 

The UVES data were taken during two distinct ob¬ 
serving runs: 13 hrs in Service Mode (ESO 79.A-0600) 
and 1.5n in Visitor Mode (ESO 80.A-0364). We used a 
combination of 390 + 564, 390 + 580 and 390 + 600 nm 
central wavelength settings appropriate to the range of 
wavelengths for the lines we were seeking. The total ex¬ 
posure time for each object was split into two or three 
equal observing blocks to minimize the effect of cosmic 
rays. The slit width was 1.2 arcsec, yielding a spectral 
resolution R = A/AA ^ 45 000. A 2 x 2 CCD binning 
was used for all observations. The observational set-ups 
are summarized in Table [2l 

The data were reduced using version 3.4.5 of the UVES 
pipeline in MIDAS. Master bias and flat images were con¬ 
structed using calibration frames taken closest in time to 
the science frames. The science frames were extracted 
with the optimal option. The blue portion of the spec¬ 
tra was checked order by order to verify that all were 
properly extracted. The spectra were then corrected to 
the vacuum heliocentric reference frame. The resulting 
spectra were combined, weighting each spectrum with its 
signal-to-noise ratio. To perform absorption line analy¬ 
sis, the spectra were normalized using cubic spline func¬ 
tions of the orders of 1-5 as the local continuum. In 
this paper, we present the UVES data for the four pairs 
that will be classified as pairs suitable for wind studies, 
hereafter wind-pair. The remainder will be presented in 
subsequent papers. We find important to mention that 
UVES and SINEONI data have their wavelength calibra¬ 
tions made in vacuum. 

2.3. Ancillary data 

Eor all of the galaxy—quasar pairs, we checked for an¬ 
cillary data and found two pairs with available HST 
observations imaging. The first one is J0448-f0950 
which has HST/W FFC2 (E555W filter) data from 
iLehnert et al.l (jl999L HST proposal ID 5393). The sec¬ 
ond one is J0839-1-1112 which has HST/WFFC2 (E702W 
filter) fron i iJgT proposal 655 7 (PI: Steidel) first pub¬ 
lished in (jKacprzak et al.l[2QlQ[ ). These HST data are 
discussed later and shown in Eigure[TJ 

3. RESULTS 

3.1. Galaxy Emission Kinematics 

In most cases, the PSE cannot be estimated from the 
data itself given the small SINEONI lEU field-of-view 


(8x8 arcsec^). Here, one advantage of using galaxy- 
quasar pairs is that the knowledge of the PSE can be 
determined from the quasar continuum in the data cube. 
This information is crucial for deriving intrinsic values of 
host galaxy parameters. Moreover, fitting a disk model 
to seeing -limited data requ ires good knowled g e of t he 
PSE (see iCresci et al.l (j2QQ9[ ) and lEpinat et al.l (j2Q12f )b 

Erom lEU data, it is customary to extract moment 
maps (e.g. flux, velocity and dispersion maps) from 
the emission hne(s) spectra. This is usually done on a 
pixel by pixel basis, as m ost algorithms treat the spaxels 
to be independen t fe.g. iFdrster Schreiber et al.l l2QQ9b ; 
Cresci et "ahl I2QQ9I: lEpinat et al.l I2QQ9I : iLaw et aH I2QQ7 . 
2QQ9f ) , a condition that requires high quality data with a 
high signal-to-noise ratio in each spaxel, in order to con¬ 
strain the width and centroid of the emission lines. Here, 
we avoid shortcomings of the traditional techniques by 
comparing the three-dimensional data cubes directly to a 
three -dimensional galaxy disk model using the Galpak^^ 
tool (jBonche et al.l [201^ . The algorithm models the 
galaxy directly in 3D (x, y. A), and the model is then con¬ 
volved with the atmospheric PSE and the instrumental 
line spread function [LSE]. The (intrinsic) model param¬ 
eters are optimized using Monte Carlo Markov Chains 
(MCMC), from which we compute the posterior distri¬ 
butions on each of the parameters. The form of the rota¬ 
tion curve v{r) is given by v{r) = Unax 2/7r arctan(r/rt), 
where rt is the turn-over radius and Unax the maximum 
rotation velocity. The algorithm has several advantages: 
(i) the dynamical center does not need to be fixed spa¬ 
tially, and (ii) the SNR required per spaxel for the cre¬ 
ation of 2D velocity maps is relaxed. In addition, since 
the actual PSE is well known from the quasar contin¬ 
uum, the returned parameters, including the galaxy po¬ 
sition angle (PA, which is defined by the angle between 
the celestial North and the galaxy major axis, anticlock¬ 
wise), inclination (i), size, maximum rotation velocity 
(Unax) B are intrinsic (or deconvolved) galaxy parame¬ 
ters. Extensive tests presented in Bouche et al. (2014, 
submitted) show that the algorithm requires data with 
a SNRmax > 3 in the brightest pixel. Eor high SNR, all 
parameters can be well recovered, but in low SNR data, 
degeneracies can appear: for instance between turn-over 
radius and Unax- 

® Since the three-dimensional disk model is inclined, the value 
Umax is the deprojected maximum velocity, corrected for inclina¬ 
tion. 
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TABLE 2 

Summary of UVES observations. 


Target 

setting Ac (nm) 

Texp (s) 

Run ID® 

Date 

J0147+1258 

390+580 

4440 

SM/079.A-0600(B) 

2007-07-23 2007-08-14 

J0226-2857 

390+580 

9000 

SM/079.A-0600(B) 

2007-07-24,27 2007-09-04 

J0302-3216 

390+564 

5430 

SM/079.A-0600(B) 

2007-08-02 

J0448+0950 

390+564 

13200 

VM/080.A-0364(A) 

2008-01-28,29 

J0822+2243 

390+564 

7200 

VM/080.A-0364(A) 

2008-01-29 

J0839+1112 

390+564 

13200 

VM/080.A-0364(A) 

2008-01-28 

J0943+1034 

390+580 

9000 

SM^79.A-0600(B) 

2007-04-18,22 

J1422-0001 

390+564 

9000 

SM/079.A-0600(B) 

2007-04-12,14 

J1441+0443 

390+600 

8100 

VM/080.A-0364(A) 

2008-01-28,29 

J2357-2736 

390+564 

4440 

SM/079.A-0600(B) 

2007-05-15 


® SM stands for Service Mode and VM for Visitor Mode. 


In order to first assess the flux profile properties, expo¬ 
nential vs. gaussian surface brightness profile and axis- 
ratios, we analyzed the collapsed cubes (i.e. line inte¬ 
grated, continuum-subtracted) with the galfit2D tool. 
This tool is our custom 3D equivalent to Galfit from 
(jPeng et al.l l2QlQf ) , which fits isophotes to the images 
(with the PSF convolution) and then uses these isophotes 
to compute the radial surface brightness profile. With 
the results from the 2D algorithm, we obtain an initial 
indication of the galaxy inclination from the axis-ratio 
and the profile shape (exponential vs. gaussian), before 
analyzing the kinematics in the 3D data. We find that 
seven galaxy surface-brightness profiles can be described 
by Exponential profiles while three are best described by 
a Gaussian profile (e.g J1422-0001). 

Using the GalPaK^^ tool, we fit the kinematics di¬ 
rectly to the data-cubes. The results are shown in Ta¬ 
ble [3l We emphasize that the surface-brightness profile 
breaks the common inclination-Wiax degeneracy in kine¬ 
matic analysis. For every galaxy, we checked that the 
MGMG chain converged for each of the parameters and 
estimated the uncertainties from the last 60 percents of 
the iterations. For J1422—0001, some of the kinematic 
parameters remain unconstrained, because the rotation 
curve appears shallow such that the turn-over radius Vt 
and the circular velocity Unax are degenerate. The pa¬ 
rameters relevant for this study for defining the kine¬ 
matic major axis (PA) are well constrained, however. 

A s we will see in section 3.4 (also illustrated in Fig¬ 
ure will . galaxy inclination is a critical parameter for 
t he w ind model (details on the model are described in 
§ HU). We cross-check the inclination measured using 
various methods (mainly galfit 2D and GalPaK^^) and 
from the SINFONI and HST data sets when present. In 
particular, for J0943+1034, the galaxy’s inclination is 
set to the value obtained from the 2D profile fitting since 
GalPaK^^ did not converged for the turn-over radius pa¬ 
rameter. For J0839-T1112, the galaxy’s inclination is set 
to the value obtained from the archival HST image. For 
J0448+0950, the galaxy’s inclination obtained from the 
SINFONI and HST data are consistent. 

In Figures [T] we present the data, ie. the observed 
flux and kinematics maps, and the fitted model. For 
each galaxy, the SINFONI data are shown in the first 
row, along with the HST image when available. In the 
second row, we present the results from the 3D kinematic 
fitting with the GalPaK^^ algorithm, where we show the 
dispersion, velocity, flux and residuals maps, from right 
to left. The residuals maps are generated from the 


residuals cubes which are just the difference between the 
data and the model normalized by the pixel noise. The 
2D maps show the mean of the residuals in each spaxel 
normalized by its standard error. All the panels have the 
north up and the east to the left. 

3.2. Redshifts 

An accurate systemic redshift is crucial to charac¬ 
terize the outflow velocity and ultimately the outflow 
mass loading factor. The GalPaK^^ algorithm outputs 
the wavelength of the Ha emission line from the axi- 
symmetric disk model. Since the galaxy distribution may 
be somewhat asymmetric, this sometimes lead to a red- 
shift bias. 

Therefore, we use two different methods. First, we 
determine the redshift from the mean of the wavelength 
of the reddest and the bluest (gaussian) Ha emission lines 
along the kinematic major axis. As a second check, we 
create a pseudo-longslit along the kinematic major axis, 
and determine Zsys from the sharp transition in the p- 
V diagram. We find that both methods yield consistent 
results. These redshifts are listed in Tabled 

The resulting intrinsic galaxy parameters will now al¬ 
low us to build a cone model in order to reproduce the 
data for galactic outflows. 

3.3. Azimuthal dependence 

In order to begin the wind modeling, we must first 
select galaxy-quasar pairs for which the quasar LOS in¬ 
tercepts the galactic winds. This can be achieved using 
the quasar azimuthal angle a between the galaxy major 
axis and the quasar (Figure [3j), because the presence of 
strong Mg II absorbers is a strong fu nction of a as demon¬ 
strated by numerous recent studies (iBordoloi et al.|[2QTl]: 
iBouche et al.l I I 2 QI 2 I: IBordoloi et al.l[^l4lLan et al.ll2Ql^ 

Hence, we use the quasar position relative to the associ¬ 
ated galaxy major axis, using the inclination and major- 
axis determined from the SINFONI data (Figs[T]), to clas¬ 
sify the different galaxy-quasar pairs in two main cate¬ 
gories (pairs suitable for wind studies (wind-pair) for 
likely outflows and pairs suitable for accretion studies 
(inflow-pair) for likely inflows). 

Figure [4] shows the azimuthal angles a versus the 
galaxy’s inclinations for our SIMPLE sample of 10 
galaxy-quasar pairs. Pairs with 60° < a < 90° are se¬ 
lected to be wind-pairs. Pairs with 0° < a < 30° corre¬ 
spond to the cases where the quasar LOS does not probe 
outflows but rather probes the extended parts of gaseous 
disks, where the gas can p otentially (or i s likely ) to be 
inflowing (inflow-pair) as in lBouch6 et al.l (j2013f i. 
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TABLE 3 

Kinematic and morphological parameters. 


Galaxy 

b (kpc) 

a (") 

inclination (°) 

PA (•) 

flux 

4^ ax 

redshift 

r 1 

Profile 

Glass 

(1) 

(2) 

(3) 

(4) 

(5) 


(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

J0147+1258 

I7.9TI.02 

30T30 

24.4T3.3 

-69T3 

1.63 

• 10-1® 

241T38 

1.0389 

7.11 ± 0.20 

EXP 

face-on 

J0226-2857 

< 2.0TI.01 

56T3.0 

47.9T1.0 

9ITI 

2.01 

■ 10-1® 

50T12 

1.0223 

2.69 ± 0.04 

EXP 

Ambig. 

J0302-32I6 

I9.7T0.95 

I6TI5 

30.4TI.5 

-37T3 

2.70 

• 10-1® 

180T15 

0.8223 

8.99 ± 0.31 

EXP 

face-on 

J0448+0950 

I3.7T0.96 

79T3.0 

52.Oil.2 

3ITI 

5.03 

• 10-1® 

253T10 

0.8391 

7.85 ± 0.07 

EXP 

wind-pair 

J0822+2243 

2I.8T0.95 

32T30 

I7.9T0.7 

I68TI 

5.04 

• 10-1® 

328T14 

0.8102 

4.14 ± 0.06 

EXP 

face-on 

J0839+III2 

26.8T0.94 

59T6.0 

72±5t 

I39T4 

1.53 

• 10-1® 

113T20 

0.7866 

5.65 ± 0.29 

EXP 

wind-pair 

J0943+I034 

24.3T1.01 

32T3.0 

43±5tt 

I40TI 

3.81 

• 10-1® 

327T10 

0.9956 

8.73 ± 0.21 

EXP 

inflow-pair 

J1422-0001 

12.7T0.98 

I7T5.0 

55T5 

8IT3 

8.93 

• 10-11' 

130±20ttt 

0.9096 

4.30 ± 0.16 

GAU 

inflow-pair 

J144I+0443tt 

10.1T1.02 

90T6.0 


87T4 

6.62 

•10-17 


1.0384 

2.99 ± 0.18 

GAU 

wind-pair 

J2357-2736 

6.7T0.95 

68T4.0 

5I.6T2.2 

I09T2 

1.29 

• 10-1® 

187T15 

0.8149 

5.53 ± 0.14 

GAU 

wind-pair 


(1) Quasar name; (2) Impact parameter; (3) Azimuthal angle a (Section 13.31) : (4) Galaxy inclination (degrees); (5) Position Angle 
(degrees); (6) Integrated Ha flux of the galaxy (ergs“^ cm“^); (7) Maximum rotation velocity (kms“^); (8) Hq; redshift (see ^ 13.21) : (9) 
Half-light radius (kpc); (10) Assumed flux profile (Exp. or Gau.); (11) Glass (inflow-pair/wind-pair) based on a selection. the 
inclination is determined from the HST data. galaxy parameters are derived from 2D fitting (galflt2D). turn-over radius is fixed to 

nlrii2 = 0.25. 


Naturally, the azimuthal angle for galaxies with low 
inclinations, corresponding to face-on cases, is very dif¬ 
ficult to constrain. These are then indexed as ’face-on’ 
cases. Pairs with a ^ 45° correspond to ambiguous cases 
where it is difficult to argue for outflows or inflows. For 
instance J0226-2857 falls into that category with the ad¬ 
ditional difficulty that this galaxy has a very low impact 
parameter {b = 0.3 arcsec or <2.0 kpc), i.e. the line-of- 
sight is likely dominated by absorption from the galaxy 
interstellar medium. 

Figure 0] shows that four galaxies are favorable to 
study galactic winds properties: J0448+0950, J2357- 
2739, J0839-M112 and J1441+0443, and are classified 
as wind—pairs in Table [3l J1441+0443 is excluded from 
subsequent analysis because our SINFONI data does not 
meet the requirement of SNR^ 3 imposed by our inten¬ 
sive tests of the GalPaK^^ algorithm. 

4. WIND PROPERTIES ANALYSIS 

4.1. Wind sub-sample analysis 

For each galaxy-quasar pair, we have the quasar 
spectrum for all the SIMPLE galaxies taken with the 
VLT/UVES instrument. In these spectra, we identified 
three main absorption features: the Mg ll (A2796 A2803) 
doublet and the Mgl (A2852) absorption line. Because 
of our selection in of 2 A (to ensure the host was 

in the SINEONI field-of-view), the Mgll doublet is sat¬ 
urated. The Mgl absorption line is not saturated and 
in most cases shows an asymmetric profile. This asym¬ 
metry can also be seen in the Mgll doublet, but less 
obviously. We center the spectrum on every absorption 
line using the derived redshifts. Eor each absorption 
line, we transform wavelength to velocity, taking the vac¬ 
uum wavelengths. Erom the absorption system kinemat¬ 
ics and geometrical properties of the galaxy, we can now 
build the winds model for the three galaxies. 

4.1.1. Cone wind-modeling 

We follow iBonche et al.l (|2Q12[) and iKacprzak et al.l 
(|2Q14D in modeling the wind as a bi-conical outflow us¬ 
ing the geometric parameters (inclination, a) set by the 
SINEONI data. The principle is to create a cone perpen¬ 
dicular to the galactic plane, fill it with uniformly dis¬ 
tributed particles and assume that the mass flux is con¬ 


served. The particles represent cold gas clouds entrained 
in the wind, since the equivalent width of the absorption 
l ines is the snm/combinatio n of several saturated lines 
(|Menard fc Chelonchel[2QQ^ . each of which corresponds 
to a ‘cold’ gas cloud {10^K) entrained by supernovae 
heated hot winds (T > 10^K). Since the galaxy inclina¬ 
tion, PA and azimuthal angle are previously determined 
from observations using GalPaK^^ and other methods, 
the only free parameters are Wut and the cone opening 
angle 6>max (see the effect of both in the Appendix). Eor 
simplicity, we assign the clouds a constant radial veloc¬ 
ity Wut, i-e. we assume that the LOS intercepts the 
clouds far from the acceleration region. 

The cone is built along the x, y, z axes: x and y repre¬ 
sent the sky plane and z corresponds to the cone height. 
Eor a galaxy with 0° inclination, the cone direction will 
be along the line of sight. We then rotate the cone along 
the y-axis to match the galaxy’s inclination derived from 
our SINEONI data and create a simulated absorption 
profile from the distribution of cloud velocities projected 
along the quasar LOS (z axis). 

We generate ^10^ particles in a cone, which are 
grouped by bins of projected velocities. The quasar 
LOS is set by the impact parameter (6) and a, both of 
which are derived from the SINEONI data cubes. Due 
to the Monte Garlo generation of particles, stochastic ef¬ 
fects create noise in the simulated profiles. This noise 
does not impact the resulting equivalent widths and thus 
the derived outflow velocities. 

We then convolve the particle velocity distribution 
with the UVES instrument resolution. Additionally, in 
order to simulate the instrument noise, we add a ran¬ 
dom Poisson noise to the simulated profile. This random 
Poisson noise has the same signal to noise ratio as the 
data and provides for a more meaningful comparison. In 
order to give an intuitive feel for this geometric model, 
we show in the Appendix examples of simulated profiles 
using different galaxy inclinations, outflow velocities and 
opening angles. 

4.1.2. Galaxy contribution model 

Since our sample consists of pairs with small impact 
parameters (6 < 20 kpc) and with inclined galaxies (from 
^ 18° to ^ 55°), we improve our model by adding the 
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Fig. 1. — From left to right: top: The HST/WFPC2 image (when available), the Ha J-band SINFONI flux (ergs“^ cm“^), the velocity 
map (in kms“^) derived from the SINFONI data and the dispersion map (in kms“^). bottom: the residuals cube represented in 2D (in 
cr), the intrinsic reconstructed galaxy with GalPaK^'^ (deconvolved from the PSF given by the quasar), its velocity map and the dispersion 
map. The quasar position is represented by the white contours on the observed flux maps when present in the map. In each panel, north 
is up and east is to the left. 







































































7 


J0448+0950 


HST F555W image 


SINFONI Ha flux 




Vel. map 


Disp. map 


1 

0 

-1 




0 ) 1 
o 

0 


-1 


-10 1 

Residuals 



-10 1 

6x{arcsec) 


-1 0 1 

Sx{arcsec 



-10 1 

Model Flux map 


le-18 


-10 1 

Model Vel. map 


-10 1 

Model Disp. map 


n 

-2 




-10 1 



200 


100 


■200 


100 


J0822+2243 


SINFONI Ha flux 



le-18 


-10 1 


Vel. map 


-10 1 


Disp. map 



■100 


le-18 


-10 1 

Model Vel. map 


-10 1 

Model Disp. map 



■lO 


■100 


■10 


-10 1 


HST F555W image 



-10 1 

Residuals 



J0839+1112 

SINFONI Ha flux , _ Vel. map 



le-18 

1 0 

fS). 

lioo 

50 

0 

^ ' 



-50 

-100 

i 


-10 1 

Model Flux map 


le-18 


-10 1 

Model Vel. map 


100 


-10 1 

Model Disp. map 


1 

0 

-1 




-10 1 



■80 


■50 


Fig. 1.— (continued) 
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Fig. 1.— (continued) 
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Fig. 2. — Same as Figure^ but for the J1441+0443 galaxy. This galaxy has the lowest SNR, below the threshold where we can trust 
the GalPaK^-^ results. Even if this galaxy is classified as wind-pair from its apparent PA, the low SNR does not allow us to build a wind 
model. 



Fig. 3.— Scheme showing the alpha angle corresponding to the 
angle between the galaxy major axis and the quasar position. 

galaxy contribution for the quasar-galaxy pairs with the 
lowest impact parameters (b < 10 kpc) such as J2357- 
2736 fSectiou r4.2.3p . The procedure is nearly the same as 
the cone model: we generate particles in a disk with an 
exponential distribution from the center to the edge. We 
take the galaxy half light radius derived with GalPaK^^ 
to estimate a realistic contribution from the disk. The 
thickness of the disk is set to be 0.15 times this radius. 
We assign the particles a constant circular velocity cor¬ 
responding to the maximum velocity of the galaxy. The 
velocity distribution of the disk is naturally strongly de¬ 
pendent on the azimuthal angle with a maximum offset 
at a = 0° and a distribution centered around 0 kms“^at 
a = 90°. 

4.2. Comments on individual wind-pairs 
4.2.1. J0U8+0950 

The galaxy near the quasar J0448+0950 has an impact 
parameter b = 13.7 kpc and an Ra ffux of 5.03 x 10“^^ 
ergs“^cm“^. Its azimuthal angle a of ^79° and incli¬ 
nation i of ^52° make it a wind-pair (Figure |4|). This 
galaxy has a SFR of ^13 Mq yr“^ (see section 4) and a 
redshift of 0.8390. 

In addition to our SINFONI data, we retrieved ancil¬ 


lary data from HST/WFPC2 (F555W filter). These HST 
data allow us to compare the morphology of the galaxy 
(see Figure [I]) with the SINFONI one. In both data sets, 
one sees that the galaxy has an asymmetric ffux distri¬ 
bution (Figure [1]) with a brighter area somewhat offset 
with respect to the kinematic center. Comparing the 
HST image and SINFONI ffux map (the quasar was sub¬ 
tracted in SINFONI Ra ffux), the PA and inclination of 
the galaxy are in good agreement. 

After determining the geometrical parameters for this 
galaxy, we can build a cone model as described in sec¬ 
tion 14.1.11 In F^ure [5l we compare the simulated profile 
for Mg I A2852[3to the observed absorption in the UVES 
data (right column of the figure). To generate this sim¬ 
ulated profile 0, we adjusted the outflow speed Wut and 
the cone opening angle while keeping the geometrical 
parameters fixed. The best values are an outflow speed 
Vout of 115± 10 kms“^ and a cone opening angle 6>niax of 
40±5°. The errors represent the maximum allowed range 
values for Wut and 6>max. 

We note that our simulated profile reproduces the 
asymmetry and equivalent width of the observed profile. 
Note, our model does not attempt to reproduce the depth 
of the profile since it is arbitrarily normalized. The ap¬ 
parent noise of the simulated profile is due to stochastic 
effects from the Monte Carlo particle distribution. 

Outflow rates and mass loading factors for e ach g alaxy 
identified as wind-pairs are detailed in section lT3l 

4.2.2. J0839-hlll2 

In our sample, the galaxy towards J0839+1112 has 
the largest impact parameter of b of 26.8 kpc. With 
an Ra flux of 1.53 x 10“^^ ergs“^cm“^, an inclination 
i of ^72° and an azimuthal angle a of ^59°, this galaxy 
also belongs to the wind-pair subsample defined in sec¬ 
tion [331 Its SFR is ^3.4 Mq yr“^ , and it has a redshift 
of ^ =0.7866. 

In Figure [T] we compare archival HST/WFPC2 
(F702W filter) images to our SINFONI Ra data. Both 
data sets show a slight asymmetry in the galaxy flux dis¬ 
tribution and a similar PA. For the galaxy inclination, we 
used galfit2D on the HST image to cross-check the results 
from the 3D fitting and found good agreement between 
them and found that the results were within ±15°. 

As for J0448±1112, we generated a simulated profile 
from the wind cone model using the geometrical param- 

^ We use Mg I A2852 since Mgii A2796 is saturated. . 

® The galaxy contribution is also considered in the simulated 
profile (the redshifted contribution in upper right of Figure O. 
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Fig. 4.— Galaxy inclinations for the SIMPLE sample as a function of the azimuthal angle a. Note there are three types of galaxies in 
this sample: the wind-pairs which have an azimuthal angle larger than 60zbl0°, the inflow-pairs with a lower than 60zbl0° and pairs 
that are ambiguous due to uncertainty on ct. It is difficult to derive the azimuthal angle for a nearly face-on galaxy. The wind-pair and 
inflow-pair classes describe the fact of having the quasar absorptions tracing outflows and inflows, respectively. 


eters from the SINFONI+HST data and adjusted the 
outflow speed Vout and the cone opening angle 6>, while 
keeping the geometrical parameters fixed. Figure [6] shows 
the simulated profile and the Mgl absorption from the 
UVES data (right column of the figure). We constrained 
an outflow speed Wut of 105±10 kms~^and a cone open¬ 
ing angle of ^max of 30±5°. The impact parameter b is 
too high to consider any contribution from the galaxy. 

4.2.3. J2357-2736 

The last individual case from our wind subsample is the 
galaxy along the J2357—2736 LOS. The host galaxy has 
the smallest impact parameter b to the quasar LOS with 
b of 6.7 kpc. This galaxy has an Ra flux of 1.29 • 10“^^ 
ergs“^ cm“^ and a SFR of ^3.3 Mq yr“^ . Its inclination 
i is ^51.6° and it is classified as wind-pair because of its 
azimuthal angle a of ^68°. 

As in the previous two cases, we generated a simu¬ 
lated UVES profile using the wind model described in 
section 14.11 Eigure [3 (bottom) shows the UVES Mg l 
A2852 absorption profile, whose asymmetry is reversed 
compared to the two other cases with a maximum op¬ 
tical depth at U ^ 0 kms~^. However, any constant 
win d spe ed model will have an outward asymmetry (Eig¬ 
ure ESI) and the data clearly shows the opposite, an 
inward asymmetry. Eor a profile with inward asymme¬ 
try, the strongest part of the absorption profile is located 
closer to the systemic velocity (e.g. bottom left of Eig¬ 
ure IZl). An outward asymmetry profile has the oppo¬ 
site behavior (e.g. Eigure [5]). This inward asymmetry is 


seen in the other non-saturated transitions (Zn ll(A2026), 
Mgl(A2026), Mnll(A2576, A2594, A2606)) present in the 
UVES data (Eig.®. 

Contrary to the other two cases, this galaxy has a very 
low impact parameter (6 ^ 6.7kpc), where the assump¬ 
tion of constant wind speed might break down. Indeed, 
the low-ionization material in momentum-driven winds 
and energy-driven winds is thought to be accelerated 
(e.g. iMurrav et al.|[2QQ5l: iSteidel et al.ll2QlQl) by the hot 
gas, by the radiation pressure, or both. 

Therefore, instead of using a constant wind speed, 
we added a generic velocity profile such as U(r) = 
Vout 2/7r arctan(r/ro) where r is the distance from the 
galaxy and tq is the characteristic turn-over radius. Eig¬ 
ure [3 shows the behavior of this model (dashed line) on 
the profile asymmetry for different values of vq , illustrat¬ 
ing that the asymmetry reverses as ro increases. The ac¬ 
celerated wind model that best describes the data is the 
one with ro = 10 kpc shown above the UVES spectrum 
in Eigure [3 Similar to J0448+0950, we also included a 
contribution from the galaxy which appears to account 
for the bluest components. 

Eor this case, we found an outflow speed Vout of 
130±10kms“^using a cone opening angle 6>niax of 45±5°. 
Note that there is no degeneracy between Vout and ro as 
the various simulated profiles shown in Eigure [3 are for 
the same outflow velocity. In other words, Vout is set by 
the reddest part of the profile, whereas ro is constrained 
by the profile shape. 

Having determined outflow velocities for the three 
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Fig. 5.— Representation of the cone model and quasar spectrum associated with the J0448+0950 galaxy. Left: the cone model seen in 
the sky plane where the y axis corresponds to the galaxy major axis and x to its minor one. The gray circle represents the inclined galaxy 
disk and the black circles illustrate the gas outflow cone. Bottom left: a side view of the cone where the z-axis corresponds to the quasar 
LOS direction with the observer to the left. Right: Normalized flux for the Mgl A2852 absorption line observed with UVES (bottom) where 
we can see an outward asymmetry, and the reconstructed profile (top). The red dashed line gives the simulated profile without taking into 
account the galaxy contribution. The black line does take into account this contribution. Note that this model does not reproduce the 
depth of the absorption line. 


wind-pairs, we now focus on deriving the ejected outflow 
rates Mout together with the mass mass loading factors 

T]. 

4.3. Outflow Rate 

For each galaxy, the equivalent width of the Mgl 
A2852 absorption lines only depends on ^max and Vout' 
since every particle has the same velocit y or accelerated 
velocities up to Vout (see Section r4.2.3p . the projected 
velocity depends on the particle position in the LOS. 
The equivalent width is due to the accumulation of 
the particles projected velocities. We tested several 
velocities and opening angles in order to fit this width. 
The profile asymmetry depends of the system geometry, 
the particle density and the outflow velocity ( 14 ut)- 
Quasar LOS that cross the wind cone can constrain 
the ejected mass r ate according to relation (1) from 
iBonche et al.l (j2Q12f ): 


^out ~ R • A^h(^) • b • Vout • — • ^max (1) 

^out ^ R Nu{b) b Lout ^max 

O.5M0yr“^ 1.5 10^^cm“^ 25kpc 200km s~^ 30° 

/i being the mean atomic weighL b the impact param¬ 
eter, 6>niax the cone opening angl^j and A^h(^) is the gas 
column density of hydrogen at the b distance. 

The only parameter that remains to be constrained is 
the column density A^h(^)- In order to determine the 
gas column density 7 Vh ^ we use the empirical relation of 
iMenard fc Chelonchd (j2009[ ) between neutral gas column 
density and Mgll equivalent width 

ATjjj = log[(3.06±0.55) X 10^^ X (V1/A2796)1.7±0.26]^ ^2) 

^ ^max is defined from the central axis, and the cone subtends 
an area E of tt • 0^^^. 
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Fig. 6.— Same as Figure [5]for the J0839+1112 galaxy where no galaxy contribution can be seen. 


This relation together with the tight correlation be¬ 
tween Mg II equivalent width and dust content (as 
determined statistica ll y fro m quasar extinction) from 
iMenard fc Chelonchd (|2QQ9f ) leads to a gas-to-dust ra¬ 
tio slightly smaller than that of the Milky Way H l 
column densities of log(A^jjj) =19.5 and above. Fur¬ 
thermore, the redshift evolution of the dust content of 
Mg II absorbers extrapolated to z = 0 shows that Mg ll- 
selected aborbers extend the local relation between vi- 
sn al extinction Ay an d the total hydrogen column A^h 
of lRohlin et (Il97l . This in turn indicates that the 
ionized gas contribution is negligible in regions with 
Hi columns above log(A^jjj) =19.5 , as also argued by 
iJenkinsI (j2QQ9f ) , and that one can use the correlation be¬ 
tween Mg II equivalent width and j as a proxy for the 
A^h gas column density. 

Given our selection criteria of > 2 A, we are 

very likely in a regime where the gas is mostly neu¬ 
tral. For our three wind-pairs sight lines, the Hi col¬ 
umn densities are: log(A^jjj) « 20.3 for J0448-t-0950, 
log(Arjjj) 20.1 for J0839+1112 and \og{N^^) ^ 19.9 
for J2357—2736. The rest equivalent widths de¬ 

termined from the UVES data and the corresponding A^h 
column densities for the wind-pair galaxies are listed in 


Table m In future work, we will be able to measure A^jj j 
directly from UV spectroscopy with HST/COS. 

Figure [9] shows the Mgll rest equivalent width 
as a function of impact parameter b for quasar- 
galaxy pairs where the quasar is aligned with the mi- 
nor axis (wind-pairs) from various literature samples 
(jKacprzak et al.l [2011allH ) and this paper. This figure 
shows that the tigh t anti-correlat i on be tween impact pa¬ 
rameter b and Wr (|Bonche et al.l (j2012f )) is confirmed at 
6 < 30 kpc. The solid line traces the fiduci al 1/6 relation 
fo^mass- conserved bi-conical outflows (see iBouche et al.l 

[20T11 . 

Table |4] lists our estimated outflow rates determined 
using Eq. [H In order to determine the error bars for 
Afout, we take the maximum error of every parameter 
used to derive it. We thereby objectively determine the 
maximum uncertainty on the ejected mass rate. We also 
note that the errors on Mout are dominated by the errors 
on 

Erom the outflow rates, we compute the mass mass 
loading factor r] by comparing it to the SFR. W e derived 
the SER from Ra u sing the iKennicnttI (jl998f ) calibra¬ 
tion, which assumes a lSalpeterl (jl955f ) Initial Mass Eunc- 
tion (IME): 
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Fig. 7.— Left column: Mgl absorption from the J2357—2736 galaxy indicate that the model must include accelerating winds., from top 
to bottom: simulated absorption profiles with ro=l,5,10 kpc, and the UVES data centered on Mgl A2852 . Notice that the asymmetry 
changes as ro increases, it goes from outward to inward asymmetry. Right column: The velocity profile corresponding to the associated 
simulated profile to the left where the turn over radius of the velocity profile (ro) varies, from top to bottom: ro=l,5,10 kpc. The red 
dashed line represents the distance between the galaxy and the quasar LOS (6/sin(Q;)/sin(z)), corrected for the inclination i. The final 
simulated profile is the one directly above the data, with ro=10 kpc. 
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Fig. 8.— Absorption lines observed with UVES for J2357-2736. 
The solid lines are a Gaussian smoothing of the data (with a sigma 
of 1 pixel) to aid the eye. The grey shaded area in the Znll(A2026) 
panel is the absorption expected from the Mgl(A2026) line (de¬ 
rived by scaling the observed optical depth of the Mgl(A2852) line, 
ignoring potential saturation effects). 


SFR{Mq yr-1) = 7.9 x 10“^^ Lna (3) 

where Lho is the Ha luminosity in ergs“^cm“^. We 
note that the SFRs for the Salpeter IMF with no extinc¬ 
tion co pection are identi cal to usi ng; a dust c orrec tion of 

1 mag (jZahid et al.ll2Q13l) with the iChabrierl (|2QQ3f ) IMF, 
as the two IMFs are offset by -0.25 dex (see Table 2 in 
iBernardi et al.l (j2QIQf )b 

The results for the three galaxies are shown in Table 01 
If we assume that galactic winds are symmetric with 
respect to the galactic plane (Figure [3]), the total ejected 
mass rate for a galaxy must be increased by a factor of 2 , 
which gives Mout ~ 9 M 0 yr“^ for J0448+0950, Mout ~ 

2 M 0 yr“^ for J2357-2736 and Mout ~ 7 M 0 yr“^ for 
J0839+III2. 


Considering the ejection velocity of the winds (115, 
105 and 130 kms“^ for J0448+0950, J0839+1112 
and J2357—2736 respectively), it is interesting to test 
whether this velocity is large enough for the gas to leave 
the galaxy halo or if it will end up falling back onto the 
galaxy. The escape velocity Ksc for an isothermal sp here 
is given by the following Eq. 0] (jVeilleux et al.l[^005f ): 


Fesc — 


1 Tin 


Rv 


(4) 



Fig. 9.— as a function of impact parameter b for galaxy- 

quasar pairs classified as wind-pairs. 

where Tdnax is the maximum rotation velocity of the 
galaxy and i?vir is the virial radius. Given that our galax¬ 
ies halos have a mass close to 10 ^^ M 0 , their virial ra¬ 
dius is approximately i?vir ~ Knax/ 10 M( 2 ;) where H{z) 
is the Hubble constant at redshift 2 ;. For the wind-pairs, 
their virial radii are 225 kpc for J0448+0950, 103 kpc 
for J0839-t-1112 and 168 kpc for J2357—2736. We give 
these results in Table 01 along with the results on mass 
loading factors r] = Mout/SFR. 

The ratio Wut/Wsc < 1 shows that the ejected gas does 
not escape from the galaxy halo and should therefore fall 
back into the galaxy. This gas contributes to the regula¬ 
tion of star formation in the galaxies. Other cases may 
be possible: for a galaxy with a low inclination (nearly 
face on) and a low impact parameter, the quasar LOS 
can easily track outflowing and inflowing materials. 

5. CONCLUSIONS AND DISCUSSIONS 

In this paper, we studied gas flows around star-forming 
galaxies using the SIMPLE sample of galaxy—quasar 
pairs (paper I). The galaxies in this sample are located 
within <3" (<20 kpc) of the background quasar sight- 
lines due to the selection of absorption with rest equiv¬ 
alent width > 2 A. Thanks to the SINEONI lEU on the 
VLT and the new algorithm GalPaK^^, we were able to 
recover the intrinsic morphological and kinematic prop¬ 
erties of the galaxies from their Ha emission (Eig. 1). 
The galaxies in our sample can be classified as wind-pairs 
or inflow-pairs according to the apparent location of the 
quasar with respect to the galaxy major-axis (Eig. 0]). 

With this classification, we focused the analysis on the 
properties of galactic winds for the sub-sample of four 
suitable wind-pairs, although one galaxy has a SNR too 
low for a robust morphological (inclination) measure¬ 
ment. The wind properties are constrained from the 
high-resolution UVES spectra of the minor-axis quasars. 
We sh ow that a simple cone model for galactic winds 
()4.1.ip can reproduce the morphology of the UVES Mg l 
absorption profiles. The wind properties can be summa¬ 
rized as follows: 

• Like other rece n t wor ks (e.g. iRnbin et al.l |2QI4 
iGhisholm et al.l I2QI4D . outflow velocities are 
smaller than the escape velocities, so the gas traced 
by low-ionization lines does not escape the galaxy 
halo. 

• At the lowest impact parameter (5^6 kpc), one 
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TABLE 4 

Results for galaxies J0448+0950, J2357-2736, J0839+1112 together with literature results. 


Galaxy 

b (kpc) 

log(NH(b)) 

Vmax 

Vout (kms 

^max 

SFR 

Vout 

Sit 

V 

Reference 

( 1 ) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

( 8 ) 

Is) 

( 10 ) 

( 11 ) 

J0448+0950 

13.7 

20.30T0.3 

253 ± 10 

115 ±10 

40 ± 5.0 

13.6±0.3 

4-6 

0.16 

0.70 

This work 

J0839+1112 

26.8 

20.10T0.3 

115 ±8 

105 ± 10 

30 ± 5.0 

3.4±0.2 

3.6 

0.43 

2.11 

This work 

J2357-2736 

6.7 

19.92T0.2 

186 ± 15 

130 ± 10 

45 ± 5.0 

3.3±0.2 

1-2 to.7 

0.24 

0.75 

This work 

J081420G1 

51.1 

19.07T0.2 

131 ± 10 

175 ± 25 

30 ± 5.0 

5.0t 

- 0.7 

0.63 

0.42 

B2012 

J091119G1 

71.2 

19.34T0.2 

231 ± 10 

500 ± 100 

30 ± 5.0 

1 . 2 t 

<-7 Q + 12.2 

- 4.5 

0.97 

12.9 

B2012 

J102847G1 

89.8 

18.60T0.2 

162 ± 10 

300 ± 25 

30 ± 5.0 

9.0t 

- 0.7 

0.95 

0.23 

B2012 

J111850G1 

25.1 

19.97T0.2 

116 ±10 

175 ± 80 

30 ± 5.0 

7.0t 

4 1 + 10.0 
- 1.1 

0.63 

1.17 

B2012 

J225036G1 

53.9 

19.54T0.2 

240 ± 10 

225 ± 50 

30 ± 5.0 

8 . 0 t 

4 2 

- 2.2 

0.40 

1.06 

B2012 

J1659+3735 

58.0 

18.89T0.15 

140 ± 10 

40-80 

40 ± 5.0 

4.6-15 

1.6-4.2 

0.12-0.27 

0.1-0.9 

K2014 


(1) Galaxy name; ( 2 ) Impact parameter (kpc); (3) Gas column density at the impact parameter (cm“^); (4) Maximum rotational velocity 
of the galaxy (kms“^); (5) Wind velocity (kms“^); ( 6 ) Gone opening angle (degrees) (7) Star Formation Rate (M 0 yr“^); ( 8 ) Ejected 
mass rate for one cone (M© yr“^); (9) Ejection velocity divided by escap e velocity; (10) Mass loading factor: ejected mass ra te divided by 
star formation rate (for both cones); (11) References: B 2012 : IBouche et al.l (| 2012 l j . K2014: IKacorzak et al.l (|2014l j. 

t SED-derived SFRs. 
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Fig. 10.— Gomparison of predicted mass loading factors from theoretical/empirical models (curves) with values derived from observations 
(dots and triangles) as a function of the maximum rotation al velocity. The result s from this work are represented by the blue circles. The 
red circle shows the mass loading factor for a z ^ 0.2 galaxy (jKacprzak et al.1l2014I L The triangles show the results for z ~ 0.2 galaxies from 
IBouche et al.l (| 2012 l) . The gray triangles show the galaxies with quasars located at >60kpc where the mass loading factor is less reliable 
due to the large travel time needed for the outflow to cross the quasar LOS (several 100 Myr) compared to th e short time scale of the Ha 
derived SFR (~ lOMyr). The upper halo mass axis is scaled on Idnax at redshift 0.8 from IMo fc Whiti (I2002I L 
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quasar-galaxy pair (J2357—2736) has an absorp¬ 
tion profile consistent with an accelerated wind 

(14231) . 

• Loading factors rj vary between ^ 0.65 for the 
two massive galaxies and ^ 2 for the galaxy with 
the lowest mass (Figure HO]). Our results indicate 
that the mass loading factors tend to be higher for 
smaller galaxies , in agreement with theoretical ex¬ 
pecta tions fe.g. Murray et al]l2Qlll: [Hopkins et al.l 

EoH). 

Figure [TO] also includes observational constraints on 
the mass loa d ing fa ctor from the z ^ 0.1 survey of 
Bonche et al.l (2012h ftrianglesi and 2 ; ^ 0.2 pair of 
Kacppak et al I(f2014f ) (square). For the iBonche et al.l 
(120121) sample, we used SED-derived SFRs (Table 4) 
given their larger impact parameters and longer travel 
times {b/Vout > 100s of Myr) compared to the time 
scale for Ha-derived SFRs (few Myr). The SED derived 
SERs are computed from the UV-to-IR photometry (us¬ 
ing Galex+SDSS+Wise surveys) with the Co de Investi¬ 
gatin g GALaxy Emission (CIGALE) software (|Noll et al.l 
I 2 OO 9 D . The pairs with the largest impact parameters 
(> 60 kpc) are shown in grey, since these mass load¬ 
ing factors can suffer strong biases due to even larger 
travel times (^ 300 Myr). Other measurements at higher 
redshift from stacked spectra of star-f orming galaxies 
(jWeiner et al.ll200^ iNewman et al.ll2012|f indicate an a v- 

erage mass loading factor of ^ 2 (p 

The different lines in Eigure liUJ represent vari- 


ous 

theoretical (jOkamoto et al. 

I 2 OIOI: IDave et al. 

2011 

: iHookins et al.l 120121: 

Puchwein & Sorinefe] 

2013 

: IVoffelsbererer et al.l 2014 \ 

Muratov et al. 120151: 


Barai ^al l I 2 OI 5 I) and empir ical models (jZahid et al.l 


2014 iPeeples fc Shankarl 120111 ) FI . In comparing obser¬ 


vations and models, it is important to bear in mind that 
some ye fiducial s c aling r el ations in the sub-grid p hysics 
(e.g. iDave et~al ] I 2 OIII: iPnchwein fc Springell 120131: 


IVogelsberger et M] l2014f ). while ot hers are from more 
complex nnmyical yproaches (e.g. JH opkins et al.ll2012l: 
Hopkins et al I 2 OI 3 I : iMnratov et al.l 120151: iBarai et al l 
20151 ) 1^4 . and thus are more directly comparable with 
observations. 

Currently our data does not allow us to discrim¬ 
inate between energy and momentum driven winds, 
but thanks to ongoing work at redshift z ^ 0.2 with 
Keck/LRIS (Bonche et ah, in prep.; Martin et ah, in 
prep.) and to the new generation lEU Mul ti Unit Spec¬ 
trosc opic Explorer (MU SE) on the VLT (|Bacon et al.l 
[3n0l : lBacon'enil[20Tl . we will be able to significantly 
increase the sample size and put tight constraints on the 
wind scaling relations. 
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APPENDIX 

UNDERSTANDING THE GEOMETRIC WIND MODEL 

In this appendix we demonstrate how varying different parameters within the cone model impacts the simulated 
absorption profile. For the cone model (Section Id.l.ip . we change three parameters in order to investigate the behavior 
of the simula ted p rofile: for different galaxy inclinations (i), cone opening angles (0), and wind outflow velocities 
(Kut)(Figure [ am]). The general trends are as follows. The different inclinations and cone opening angles change the 
left portion of the simulated absorption profile, while the wind velocity extends th e rig ht portion. Note that, except for 
the case of 90° galaxy incli nation , all the simulated absorption profiles in Figure [AT] present an outward asymmetry. 

We also present in Figure lA^ the UVES Mgl A2852 absorption lines of our SIMPLE sample galaxies. 


2=30° 2=60° 2=90° 



Fig. a.I. — Examples of simulated absorption profiles with different galaxy inclinations (i), opening angle (0) and wind velocities (Uout): 
while each of the simulated profiles has the same number of particles, the apparent depth decreases as each parameter increases due to 
larger velocity projections for i and 0, and larger range of velocities for Uout- top row: absorption profiles for galaxies inclined at 30, 60 and 
90 degrees with Uout= 100 kms“^and 0= 30°. The noise effect is due to the Monte Garlo distribution of particles, middle row: absorption 
profiles for wind cones with opening angles of 30, 40 and 45 degrees with V5ut= 100 kms“^and 2=45°. bottom row: absorption profiles 
with wind velocities of 50, 100 and 150 kms“^with i= 45 degrees and 6= 30°. Each simulated profile has the same amount of particles 
but show a larger velocity range due to the increasing gas speed, hence the varying apparent depths. 
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Fig. a. 2 . — Here we present the UVES Mgl A2852 absorption lines for 8 out of 10 galaxies. We do not have these absorption lines for 
J0147+1258 and J0226—2857 as they fall in the gap of UVES data. The top row corresponds to UVES Mgl centered data of the four 
wind-pairs. The bottom row presents the UVES Mgl centered data of the two inflow-pairs and the two ambiguous cases from left to right. 













































